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disk and reminiscent of ROSAT observations of the NLS1, WPVS 007. We show that this is likely 
to be a "true" Sy II, without a BLR, and speculate that in its new high luminosity state we may be 
able to witness a BLR in formation. 
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1. Introduction 

The variability of the radiation output of AGN is well established and is one of their defining 
characteristics. In single observations of Seyfert I galaxies, fluctuations in the X-ray output on 
timescales of minutes to hours is the norm (e.g. NGC 4051 [[[]] and MRK 766 [Q]). The power 
spectral density function (PSD) of these variations can be described by red noise; a smooth spec- 
trum with increasing amplitude at lower frequencies (^) . The PSD of AGN has been shown to 
have a break, where the slope flattens, occurring on a timescale that scales with the luminosity 
and by implication the black hole mass. Lower mass (luminosity) AGN have a break in the PSD 
of ~ l(T 2 days (e.g. NGC 4051 [@|) whereas the most massive AGN have breaks at ~ 100 days 
(e.g. NGC 3516 [§). This means that on short timescales, low mass AGN exhibit stronger X- 
ray variability than their more massive counterparts whereas on timescales of years the variability 
should be consistent. An additional dependence on the fraction of the Eddington accretion rate, m 
has been observed, leading to the refinement that the variability can be explained by characteristic 
pertubations in a steady-state accretion disk, whose inner radius scales with Mg#/m [||]. 

The class of objects known as NLS1, show stronger variability than Sy I and exhibit more 
rapid variations. This is consistent with the belief that they are low mass systems accreting at a 
high fraction of the Eddington rate. 

Very large changes in flux due to variations in line-of-sight absorbing material have been ob- 
served in several sources. In the Seyfert 1.8 galaxy, NGC1365, factor 10 variations, accompanied 
by strong spectral changes have been seen several times (ftjh, ^) and can be modelled by move- 
ments of individual clouds in the broad-line region. 

Finally, Bl Lacs or Blazars have been famous for their volatility at all wavelengths for a long 
time (e.g. [§). 

While there has been a lot of work on short-timescale variations in specific sources and longer 
baseline studies of bright sources with RXTE, to date there have been few systematic studies of 
the variability of large samples of objects. This is mainly due to the lack of large-area sensitive 
sky surveys. The ROSAT All Sky Survey (RASS [JT^]) is an excellent resource for AGN with 
Fq.2-2 > 3 x 10~ 13 erg s _1 cm but it can only be systematically compared with data taken from 
ROSAT pointed observations taken months or a few years after the RASS completed in 1990. With 
the advent of the XMM-Newton slew survey (XMMSL1 [|P1|]) containing observations from 40% of 
the sky with fb.2-2 > 6 x 10~ 13 erg s _1 cm~ 2 , a sensible comparison of AGN flux over a baseline 
of 3-19 years can now be achieved over 40% of the sky. In this paper we make a systematic 
comparison of the long-term variability of all classes of AGN detected either in the ROSAT or 
XMMSL1 surveys. A ACDM cosmology with (co M , co x ) = (0.3,0.7) and H =70 km^sec" 1 Mpc" 1 
has been assumed throughout. 

2. Method of ROSAT/XMM comparison 

We have made a sample of AGN which have been observed (but not necessarily detected) in 
both the RASS, or ROSAT pointed observations, and XMMSL1. We define three groups of objects: 

• XMMSL1 sources, identified as AGN or galaxies, which have counterparts in ROSAT. 
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• XMMSL1 sources, identified as AGN or galaxies, which have been observed but not detected 
by ROSAT. 

• RASS sources, identified as AGN or galaxies, which have been observed but not detected in 
XMMSL1. 



We have used sources from the XMMSLl-delta3 CLEAN catalogue, released in August 2009, 
containing slew data taken between August 2001 and January 2009. 

Sources which are detected with extended X-ray emission in either XMM-Newton or ROSAT 
have been excluded. 



2.1 Upper limits 

For the XMMSL1 AGN without a ROSAT counterpart a 2-sigma upper limit to the ROSAT 
count rate has been found using the EXSAS software package ([12]). 

For RASS AGN, observed but not detected in XMMSL1, we have calculated a 2-sigma upper 
limit to the XMM-Newton count rate using a web-based upper limit server 1 ([13]). 

The final sample contains 1038 AGN of which 689 were detected in both surveys, 223 have a 
RASS upper limit and 126 have an XMMSL1 upper limit. 



2.2 Count rate to flux conversion 

All XMM-Newton slew observations have been made with the EPIC-pn detector using the 
Medium filter. In Figure |l| we plot the XMMSL1, 0.2-2 keV count rate against the RASS count 
rates. We find the mean ratio of the XMM to RASS count rates is 6.94+/-0.20. From now on we use 
a factor 7 to normalise the count rate ratio, such that an intrinsic ratio of 7 is considered to show the 
same 0.2-2 keV flux in the two missions. The conversion factor from count rate to flux is crucial 
for comparing fluxes from different insuuments and we need to be very careful when considering 
the effect of the source spectrum on this conversion. In Figure ^| we show the count rate ratio 
(normalised by the factor 7; Rj) as a function of the source spectrum. In principle there is a large 
variation of Rj if the absorbing column is very low or very high. With these AGN observations, 
NH is practically limited to a minimum of ~ 1 x 10 20 cm~ 2 , due to our Galaxy, and a maximum 
of ~ 10 22 cm~ 2 due to the flux limitations of the two surveys 2 . This can be seen clearly in the 
sky distribution of the objects in the sample where few objects lie in the Galactic plane (Fig. 
Similarly we can infer that the intrinsic absorption can not much exceed ~ 10 22 cm~ 2 . In this range, 
and for spectral slopes of 1.0 to 3.0, we see that a source with equal flux in RASS and XMMSL1 
observations will have 0.3 < Rj <1.5. 

We now define the variability ratio, Ry, to be the ratio between the higher and lower count rate 
irrespective of the observatory 3 

'http://xmm.esac.esa.int/external/xmm_products/slew_survey/upper_limit/uls.shtml 

2 This applies to cold absorption. The behaviour of ionized absorbers which permit the passage of low-energy X-rays 
is quite different 

3 Hence Ry = 1 when R-; = 1 and Ry = 2 when Rj = 2 or Rj =0.5. 
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Figure 1: A comparison of the ROSAT and XMM-Newton (0.2-2 keV) slew count rates 
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Figure 4: The X-ray variability of our sample as a function of XMMSL1 0.2-2 keV luminosity; left panel: 
individual sources, right panel: median variability per luminosity bin. 



3. Results of ROSAT/XMM comparison 

From the full population of 1038 AGN we find a median count rate ratio (Ry) of 2.11, with 
68.4% of the sources having R v < 3 and 4.9% with R v > 10. 

We have made simulations of a non-varying population of sources with the Galactic Nh ran- 
domly selected from the actual sky values of the 1038 AGN, power-law indices selected from the 



population of ROSAT AGN slopes reported in Q14| ] (T = 2.05 ±0.55) and count rates randomised 



using the actual errors in the XMMSL1 and RASS measurements. The simulated sources have a 
variability ratio with a median value Sy=l. 295 ±0.003. 

Following, [16] and [|l7|] we calculate the true variability signal (V) by: 



V 2 = (R v -l) 2 -(S v -\) 2 (3.1) 

where Sy is effectively the noise due to the spectral effects and measurement errors. This 
gives a median variability, V, of 107 ±7% in the 0.2-2 keV energy band. In the UV band the mean 
long-term variability was found to be 35% from a survey of 9000 SDSS quasars ((HHl). This result 
then continues the trend of higher photon frequencies exhibiting larger amplitude variability, seen 



on shorter timescales (e.g. NGC 5548; [|19J). 



3.1 Luminosity dependence 

In Fig. ^ we plot Ry against the X-ray luminosity for the sources where the redshift is known. 
No significant correlation is found. As we are probing variability timescales far longer than the 
PSD break timescale it is normal that the luminosity-dependence, seen at shorter timescale, has 
no influence. This result shows that no further luminosity (mass) dependent effect is affecting the 
longer-term variability. 

3.2 Source category 



We have cross-correlated our sample with the Veron catalog of Quasars & AGN [22] to find 



the AGN class for each source, where known. In Table 1 we separate the sources by AGN class and 
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Table 1: Variability statistics for each AGN class 



CLASS" 


No. b 


Median 6 ' 




Fraction (%) d 






Var. (%) 


R v < 2 


R v < 3 


Rv > 10 


ALL 


1038 


107 ±7 


46.8 


68.4 


4.9±0.8 


SI 


318 


90 ±12 


51.8 


72.0 


4.1 ±1.4 


S1.2 


48 


66 ±27 


62.5 


78.9 


0.0 


S1.5 


83 


112±26 


45.1 


65.7 


3.3±2.4 


S 1.8/1.9/2 


37 


155±44 


34.3 


55.9 


13.3±6.0 


QSO 


232 


130±18 


42.4 


62.1 


4.8±2.0 


NLS1 


64 


99 ±28 


48.1 


72.5 


4.1 ±2.9 


Blazar 


142 


109 ±19 


43.9 


73.4 


3.4±1.7 



" AGN class: Seyfert 1, 1.2, 1.5, a combination of 1.8, 1.9 and 2, QSO, NLS1 and a combination 
of Blazars, blazar candidates and highly polarised quasars. 
h The total number of sources in this category. 

c Median variability, V, expressed as a percentage, after correcting for spectral effects (see text). 
d The percentage of sources with Rv less than 2, less than 3 or greater than 10. 



calculate the median variability and the fraction of sources showing Ry < 2, < 3 and > 10. The 
calculations ignore sources with upper limits that fall outside of these limits. 

Seyfert II galaxies have the highest percentage of sources showing strong (Ry > 10) variability. 
It is not clear whether this is due to changes in line-of-sight absorption or changes in the intrinsic 
emission and these data can not help distinguish between the two possibilities. A closer look at one 
particular example is needed. 



4. GSN069 

During a slew on July 14, 2010 XMM-Newton detected soft X-ray emission at 1 .5 ±0.3 counts 
s _1 from a position consistent with the nucleus of GSN 069 (also known as 6dFg01 19087-341 131; 
z=0.01816) (see Fig. g). Previous ROSAT survey and pointed observations failed to detect the 
source and give upper limits, to the unabsorbed flux, 30-360 times fainter than the XMM-Newton 
slew detection (Table ||). 

Two, consistent, optical spectra, taken in 2001 and 2003 in the 2dF and 6dF surveys [ j20| ] show 
unresolved Balmer lines (FWHM < 200 km s _1 ) and line ratios consistent with a Sy 2 classification 
(Fig. |). 

The source has been monitored with the SWIFT-XRT since its discovery. The SWIFT obser- 



vations have been analysed following the procedure outlined in Evans et al. (2009) [21] and show 
stable soft X-ray emission to within a factor 2 (Fig. |7|). 
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Figure 5: A DSS image and contours of GSN 069: the red circle represents the XMM-Newton slew error 
circle, blue is for the SWIFT-XRT position and green is the SWIFT position enhanced using UVOT 
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Figure 6: Optical spectrum of GSN 069 from the 6dF survey 
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Figure 7: The X-ray light curve of GSN 069 
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Table 2: GSN 069 observation log 



Mission 13 


Date 


Count rate 


Mux 


RASS 


1990 


< 0.0099 


< 0.11 


ROSAl-PSPC 


1993-07-13104:38:55 


- / \ f\f\ inn 

< 0.00177 


< 0.019 


ROSAT-PSPC 


l994-06-29T05:59:5l 


< 0.00129 


< 0.014 


YA/TA/r elf,,,, 

Aiviivi siew 


omn m iatho-as-ia 
zu lu-u / - m iuu.4o. i^t 




J.l ± 1.0 


SWIFT 


2010-08-27T06:05:40 


0.041 ±0.005 


4.3 ±0.6 


SWIFT 


2010-10-27T05:32:30 


0.037 ±0.004 


3.9±0.6 


SWIFT 


2010-11-24T01:19:48 


0.033 ±0.004 


3.5±0.5 


XMM pointed 


2010-12-02T10:44:18 


1.03 ±0.01 


4.32 ±0.04 


SWIFT 


2010-12-22T15:05:03 


0.029 ±0.004 


3.0±0.5 


SWIFT 


2011-01-19T06:15:29 


0.028 ±0.004 


2.9 ±0.4 


SWIFT 


2011-02-16T05:33:31 


0.029 ±0.003 


3.0 ±0.4 


SWIFT 


2011-04-25T05:00:28 


0.027 ±0.003 


2.8±0.4 


SWIFT 


2011-05-23T04:06:01 


0.035 ±0.004 


3.7±0.5 



" XMM-Newton, EPIC-pn camera: slew observation performed in FullFrame mode with the 
Medium filter; pointed observation performed in FullFrame mode with the Thin filter. SWIFT-XRT 
in pc mode. 

b counts s _1 in the band 0.2-2 keV. 

c Unabsorbed flux, Fo.2-2 keV, units of 10~ 12 erg s _1 cm~ 2 , calculated from a black body model 
with kT=58 eV and Galactic absorption (2.48 x 10 20 cm~ 2 ) 
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Figure 8: Residuals to a diskbb model fit to the XMM-Newton EPIC-pn (black), MOS-1 (red) and MOS-2 
(green) data 
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4.1 X-ray spectral fits 

An XMM-Newton TOO, with nominal exposure lOks, was performed in December 2010, 
using all of the X-ray detectors and the optical monitor. The data were reduced using XMM- 



SAS vll.O (Q15|]). EPIC source spectra were extracted from a circle of 30 arcsecond radius and 
background taken from a nearby source-free region. The two RGS cameras were combined into 
one spectrum. All the camera spectra were grouped so as to contain a minimum of 25 counts 
per bin and to oversample the instrument resolution by a maximum of 3. Fits to the individual 
instruments gave consistent results and so all spectra were fit simultaneously, with a constant used 
to compensate for the small normalisation differences. 

An initial fit with a power-law, absorbed by the Galactic column (2.48 x 10 20 ; LAB Map - 



[g3J) gave a poor fit but showed that the overall spectrum is very soft (T = 6.2). To model the shape 
better we fit the soft X-rays with a black-body or multi-colour disk model (diskbb in XSPEC) and 
the weak emission at > 2 keV with a power-law (with T fixed to 1.7). The black-body and diskbb 
models yielded better fits with = 1 .4, 2.0 respectively (Tab. ||). To fit the remaining residuals we 
tried various absorption models. Cold absorption, in the rest frame of the source, is excluded with 
an upper limit of 9 x 10 19 cm -2 . Good fits were obtained with a simple edge of energy 650 eV or 
with an ionized absorber with /V//=4-l 1 x 10 23 cm~ 2 , t, = 100 — 800 and a covering fraction of 74- 
100% redshifted by either 0.15 or 0.2 in the rest frame of the source, depending on the continuum 
model (Fig. ||). 

The unabsorbed soft X-ray luminosity is Lo.2-2 = 5.3 ±0.1 x 10 42 erg s _1 for the best fit 
model (diskbb with edge). This model greatly underpredicts the flux observed in the OM UVW1, 
UVM2 and B filters. However, the optical images of the galaxy appear extended in all filters and 
so the count rates are clearly dominated by emission from the stellar population. We therefore 
extrapolate the soft X-ray spectral model over optical, UV and EUV frequencies to obtain Lb i = 
1.8 ±0.1 x 10 erg s . The emission above 2 keV is very weak, 1.1 ±0.9 x 10 counts s _1 in 
the EPIC-pn detector, and the conversion of this count rate to an unabsorbed flux is highly model 
dependent. We find a highest, 90% confidence, upper limit of L2-10 < 8.5 x 10 40 erg s _1 , for the 
black-body with partial covering model. 

A spectral fit of the combined SWIFT-XRT data gives very similar spectral parameters to the 
XMM-Newton spectrum showing that the soft spectral shape of GSN 069 is long lasting. 

A possible explanation for the steep X-ray spectra of GSN 069 is provided by an ionized 
absorber which allows very soft X-rays to pass while strongly blocking higher energies. This model 



was used to explain the very steep spectra of WPVS 007 seen during the ROSAT era y24fl. To test 
for this we have attempted to fit the GSN 069 XMM-Newton spectrum with intrinsic power-law 
emission absorbed by up to three ionized, partial covering, absorbers. No good fit was obtained 
with a power-law model with T < 4 and therefore we conclude that the very steep soft X-ray 
spectrum is thermal in origin. 



4.2 Discussion 

The lack of a cold absorber in the X-ray spectrum of GSN 069 is difficult to reconcile with 
the classical interpretation of a Sy II where the broad lines are removed by absorption. From the 
6dF optical spectrum we find a ratio of [OIIIJ5007 / Hp ~ 1 1 which is symptomatic of a missing 
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Table 3: GSN 069 spectral fits to an XMM-Newton TOO 



Low-energy model Intrinsic absorption ^ 2 /dof 



Plaw Bbody Diskbb Edge" zxipcf* 



T kT(ev) kT(eV) Norm 


£(eV) Tau Nh % cf(%) z 




6.2 

56.3^ 
- 58.2+^ - 
62 2+ 07 

67.1 43695 
71.71™ 28184;™ 
77 4 +i.i 21 105+ 340 

"■^-2.0 ziiuJ -460 


651+^2 o.80 +Q Q f 9 - 

- 98ii« 2.7+0+ g2 +4 167 +0.005 

650+1? - 

9 8+ 01 Q9+ 8 H 
J -10 z -°-0.6 ^ -8 u ' ZZJ -0.007 


2015/107 

152/107 

116/105 

114/103 

215/107 

112/105 

116/103 



All fits included absorption by the Galactic column (Nh = 2.48 x 10 20 cm~ 2 ) and a power-law with 
slope fixed at T = 1 .7 to model the very weak emission beyond 2 keV. Errors are 90% confidence. 
a The edge model in XSPEC. b A partially ionized, partial covering model (zxipcf in XSPEC) with 
parameters, Nh in units of 10 22 cm~ 2 , log of the ionization parameter, covering fraction and redshift. 



BLR rather than absorption in Sy II [31]. From this we conclude that the BLR is absent rather 
than obscured in this case. In other words, GSN 069 appears to be another example of an AGN 
classified as a Sy II galaxy but with little or no intrinsic absorption in the X-rays and with optical 
line ratios inconsistent with heavy extinction. The X-ray spectrum could be reconciled with a Sy 
II classification if the AGN was Compton-thick in the X-rays (although the optical line ratios do 
not support this idea). However the soft X-ray variability on both short (doubling time of 800s in 
the XMM-Newton observation) and long (see Table 2) timescales rules out this hypothesis, making 
GSN 069 a "true Sy II" candidate, i.e. a candidate AGN with no BLR. 

In the paradigm where the BLR is created from outflows of material from the disk, the outflow 
may not be sustainable and hence the BLR may be absent if Lboi < 10 42 erg s _1 [26] or if the 
accretion rate is below a critical value m < 0.01 [25]. 

The black hole mass, Mbh, can be estimated from the normalisation of the diskbb model if 
the emission is thermal, as it seems to be. Following Yuan et al. 2010 [|3(J, and taking the disk 
inclination as i=40° for simplicity, we derive the black hole mass as Mbh = 1.7 x 10 5 M for a 
non-rotating black hole and Mbh = 8-4 x 10 5 M for a maximally spinning Kerr black hole. The 
accretion rate may be derived from the z-corrected temperature, to be m = 0.017 — 0.38 depending 
on the black hole spin. An independent estimate of Mbh may be found from its relationship with 
the bulge K band luminosity ft2$j. The 2MASS extended catalogue gives m^=12.75, for the whole 
galaxy, which implies Mbh < 2.5 x 10 M . The low mass of the black hole is also supported, to 
some extent, by fast variability (doubling time of 800s) seen in the XMM-Newton light curve. 

From the relationship between Mbh, m and the FWHM of the BLR lines ^ we would expect 
to be seeing broad lines with FWHM=560-1800 km s~ ! in the optical spectrum. 

The luminosity of the [OIII]5007 line in the 6dF spectrum is 1.1 x 10 40 erg s _1 which implies 



10 



Long-term AGN variability 



Richard Saxton 



a historical L/ x ,/ ~ 10 42 erg s _1 using a bolometric correction factor of 88-142 [29]; an order of 
magnitude below the current value. If radiating at an efficiency v\ =0.1, the historical accretion rate 
will have been m = 0.009 — 0.047. In the past then, the source has been accreting and radiating close 
to the critical values necessary to produce outflows and plausibly has been too weak to produce the 
BLR. With the current AGN parameters an outflow should be sustainable, in which case we can 
estimate the time necessary to form the BLR as t = v/d, where v is the velocity of expanding 
material and d is the distance of the BLR which will be ~ 10 days for Ltoi = 10 43 erg s _1 
Hence for outflow velocities of 1000-6000 km s _1 (e.g. [27]) BLR material will start to build up 
after 2-12 months. A regular monitoring of the optical spectrum of GSN 069 may be able to detect 
this process. 

The apparently redshifted, strong edge observed at ~ 650 eV may be produced by a failed disk 
outflow or from an inflowing accretion disk overlapping our line of sight. 

In summary, the soft emission from this SY II appears to be thermal in origin and its variability 
appears to be intrinsic to the central engine and not related to changes in line-of-sight absorption. 
This mechanism may also be the cause of the strong variability seen in the other Sy II galaxies, 
three of which have Lx < 10 43 erg s _1 in their high state. 
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